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We discuss the universal critical behavior of a selfinteracting scalar field theory at finite 
^ . temperature as obtained from approximate solutions to nonperturbative renormalization 

^ I group (RG) equations. We employ a formulation of the RG-equations in real-time formalism 

which is particularly well suited for a discussion of the thermal behavior of theories which 
are weakly coupled at T = 0. We obtain the equation of state and critical exponents of the 
theory with a few percent accuracy even for a relatively simple approximation to the exact 
renormalization group equations. 
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1. Field theory at nonvanishing temperature has in recent years attracted a lot of interest. It 
is relevant on widely separated scales, starting from the physics of the early universe, where 
phase-transitions may play a role in inflation, the formation of density inhomogeneities and the 
creation of the matter-antimatter asymmetry observed today, via the by now also experimentally 
studied physics of the quark-gluon plasma to critical phenomena in solid state systems. Thermal 
energies involved range from 10^^ GeV down to a few meV. On the other hand, the number of 
open questions is still large compared to other fields of theoretical physics. Even for systems that 
can be treated as equilibrium systems, there are a number of technical problems that have to be 
overcome. This holds already for very simple models like the selfinteracting scalar theory near its 
critical temperature studied in this letter. 

Perturbation theory at nonvanishing temperatures is plagued by infrared problems. This is 
mainly due to the fact that the dynamics of the theory at high temperatures is dominated by 
classical (thermal) fluctuations, which may be described by a three-dimensional effective theory. 
The infrared behavior of the three-dimensional theories prevents us from studying the critical 
behavior in straightforward perturbation theory. Indeed, simple one-loop resummed perturbation 
theory predicts a first order phase-transition for the 0(A^)-symmetric scalar theory [0] - a result 
that is known to be wrong. 

There are a number of approaches available that cure the problems of straightforward per- 
turbation theory. Most of these methods use variants of the renormalization group (RG) idea. 
The "environmentally friendly RG" (0 and references therein) is closely related to the field the- 
oretic (Callen-Szymanzik type) RG and as such in essence perturbative, avoiding however the 
above-mentioned infrared problems. This makes the method very well suited to study the critical 
behavior of systems where the basic degrees of freedom do not change dramatically through the 
phase transition, as demonstrated e.g. by the results for the critical exponents of scalar theories 
like the one studied here 0. 

Another approach is the so-called "auxiliary mass method" suggested in and applied to the 
critical behavior of scalar theories in |^ . This method is in principle nonperturbative and thus, like 
the thermal renormalization group approach employed in this letter and discussed below, suffers 
from the need to make approximations that may in general only be tested a posteriori. Results 
for the critical indices of selfinteracting scalar theories are reported in 0] and are comparable to 
the results of the present study (see below). 

The standard approach of statistical mechanics is the Wilsonian renormalization group 0. The 
basic idea of the Wilsonian RG may be implemented in different ways. It has been introduced 
in a field-theory context by Polchinski and others mainly to provide a physical understanding 
of renormalization. In recent years this approach has been applied to different problems and is 
now also known as the "Exact renormalization group (ERG)" -approach. It takes no recourse to 
perturbation theory and has been used to study strongly interacting problems both at vanishing 
and at nonvanishing temperature. Applications to critical phenomena are commonly done - in the 
spirit of the above mentioned effective theories - in the framework of three-dimensional theories 
|0]-[|13], although some work using a formulation of thermal field theory in the "imaginary time"- 



formalism has also been done [|T4| 



Recently, an implementation of the Wilsonian renormalization group for thermal field theory in 
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the real-time formalism has been proposed [|15l (we will henceforth denote this formulation as the 
"thermal renormalization group" or TRG). Even though on a formal level the connection between 
the imaginary-time and the real-time formulation of thermal field theory is well understood |16| 
and all physical quantities may be computed in any of the two approaches, if one is interested 
in time-dependent quantities it is often convenient to be able to directly work in a real-time 
approach instead of having to do analytical continuation of Green-functions which may be involved 
or even impossible if the results are only numerical (as is typically the case in applications of the 
RG). Furthermore, in the real-time formulation of thermal field theory there is a clear separation 
between "thermal" and "quantum" fluctuations. The two-point functions may be separated into 
a part identical to the propagator at vanishing temperature and a part which contributes only 
at T > 0. In the formulation proposed in |15|, only the thermal fluctuations are treated by 
renormalization group methods. The quantum fluctuations are assumed to be integrated out, i.e. 
one starts from the full physical effective action of the theory at vanishing temperature. This 
feature may seem to be a shortcoming, since it essentially prevents us from studying theories 
which are strongly interacting already at vanishing temperature^. On the other hand, having an 
exact renormalization group equation with the physical theory as boundary condition turns out to 
be a major advantage if one is interested in the thermal behavior of weakly interacting models in 
the vicinity of a phase transition which is second or weakly first order and thus not accessible in 
perturbation theory. In this case, one obtains the thermal behavior directly in terms of physical 
quantities (masses, couplings etc) instead of connecting it to the scheme-dependent parameters 
of the action of the theory. The usual formulation of the ERG defines its own renormalization 
scheme and the connection between physical couplings at T = and the parameters of the action 
has to be obtained through explicit calculation This typically involves ambiguities which are 
only fixed at the two-loop level in a perturbative language. Formal two-loop exact calculations in 
a nonperturbative scheme like the exact renormalization group are however notoriously difficult 

Another advantage of the TRG is its explicit gauge-invariance []TP[ . In the Wilsonian renormal- 
ization group approach one introduces an external scale into the theory which is used to separate 
"hard" and "soft" modes in order to be able to consecutively integrate out "momentum shells". 
This procedure violates gauge-invarianc^, and one obtains generalized Slavnov- Taylor identities 
imposing relations between the terms in the effective action Although from a formal point 
of view everything is well understood, it is often not trivial to implement these constraints in 
actual calculations. In the TRG-approach, gauge-invariance is maintained since the separation 
between hard and soft modes only concerns on-shell fields, namely the thermal fiuctuations. This 
feature is very advantageous for practical calculations, since it severely restricts the possible form 
of solutions to the renormalization group equations. 

These points make the TRG the ideal tool to study the critical behavior of (at vanishing 
temperature) weakly interacting models beyond the universal aspects, which may be calculated in 

^Although in principle the effective action of such models could be obtained by nonperturbative methods (lattice, 
ERG-equations at T = 0, etc), this is in practice a formidable task for realistic models like QCD. 

■^One may work in the background-field formulation, keeping manifest gauge-invariance with respect to the 
background-gauge transformations pffl . 
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the framework of three-dimensional effective theories. There are a number of interesting problems 
which are weakly interacting at vanishing temperature but are nevertheless not accessible through 
perturbation theory at temperatures close to the critical one. As stated above, the simplest 
examples are 0(A^)-symmetric selfinteracting scalar theories which are known to exhibit second 
order phase-transitions. 

Another important example is the electroweak sector of the standard model, which is related 
to the observed asymmetry between matter and antimatter in the universe (see |^2[ for a review). 



This model is weakly interacting at zero temperature for not too large Higgs-masses, but its 
critical behavior is not described by perturbation theory (unless the Higgs-mass is very small). 
Around the critical temperature, the model is effectively strongly interacting for a large range of 



Higgs-masses |23 



Finally, the TRG may be used to study real-time quantities like thermal damping rates |24 
etc. One could even hope to be able to access non-equilibrium quantities along these lines. 

A prerequisite for the study of thermal phenomena in the framework of the TRG is the existence 
of useful approximation schemes for the solution of the exact renormalization group equation. This 
letter presents a study of the one-component selfinteracting scalar theory as a testing-ground for 
such nonperturbative approximations. The critical behavior of this theory has been excessively 
studied with the help of the three-dimensional model. The theory is in the universality class of the 
Ising model and the critical exponents are known to high accuracy from lattice calculations, the 
e-expansion, ERG-studies etc (for a review, see |^). A quantity which is more difficult to access is 
the critical equation of state (EOS) of the theory. Recently, also the EOS has been obtained from 



lattice studies and the ERG in three dimensions (in a slightly different context, the EOS 
has also been calculated using environmentally friendly renormalization |^). We aim at deriving 
the critical behavior not from the effective three-dimensional model, but directly from the finite- 
temperature theory. The results will show that one can give useful approximations to the exact 
renormalization group equations in the case studied and will give a nonperturbative example of 
dimensional reduction. We will not be able to improve on other methods for the critical exponents 
within the approximations performed here. Nevertheless, our results are accurate to a few percent 
and do not rely on universality, thus proving the usefulness of the TRG for the study of critical 
systems and opening the way to obtain nonuniversal quantities like critical amplitudes directly in 
terms of physical quantities. 

2. Before presenting the results, we will briefly describe the formulation of the thermal renormal- 



ization group as proposed in |jT5|. The formulation of the TRG-equations relies on the closed-time 
path (CTP) approach to thermal field theory [0. In the operator formalism, one starts with 
the usual expansion of fields in terms of creation and annihilation operators a|, and a^. Operator 
averages are now thermal averages, and we have 

<a\ak'>T = {2Tc)'^2uJkN{uJk)6{k - k') 

<akai,>T = i27rf2LUk[l + Niu;k)]5ik-k') (1) 



where Uk = \Jk'^ + m? and Niojk) = {e^^/'^ ~ ^ ' order to divide the thermal modes in to 
"hard" and "soft" ones, one modifies the thermal part of these averages, i.e. the part proportional 
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to the thermal distribution function N{uJk) by introducing a cutoff-function 

{2TTf2iUkN{uJk)Qi\k\,A)6{k~k') 
{27Tf2uJk\l + N{uJk)e{\k\, A)] 6{k 



A) as 



<alak'>T 
<akCili>T 



k' 



(2) 



where G)(|A;|, A) is 1 for \k\ > A and rapidly vanishes for \k\ < A. This has the effect of suppressing 
thermal modes with three-momenta small compared to the external scale A - the external scale 
acts as infrared cutoff for the thermal fluctuations. The formulation of a generating functional for 
IPl Green-functions proceeds as usual, except that the two-point function appearing in the action 
now depends on the cutoff A. We thus obtain an effective action Fa which interpolates between 
the physical effective action r(T = 0) of the zero-temperature theory as A ^ oo (no thermal 
fluctuations included) and the full effective action of the theory at finite temperature F as A — >■ 
(all thermal fluctuations integrated out). The explicit form of the two-point function D]^^ and the 
generating functional Fa may be found in |jl5| . 

In the end of the calculation one will usually be interested in the limit A — » 0, and in the 
spirit of the renormalization group one now proceeds to deriving a differential equation for the 
dependence of Fa on the external scale. Quite analogously to the exact renormalization group 



equation proposed in [53], the thermal renormalization group equation takes the form |T5 



A|^r.|»l 




(3) 



(^ is an exact equation and one faces the usual problem of not being able to solve it without 
approximations. In the present study, we will work in lowest order of an expansion of the effective 
action in powers of derivatives. The derivative expansion effectively may be seen as an expansion 
in the anomalous dimension of the fields, and we expect it to work rather well in the scalar theory 
under study here We thus proceed by defining the functional T/^[(f)] through |^ ($i and $2 
are the physical and the "thermal ghost" field resp.) 



5F/ 



5Fa[$] 



50 

For this functional we make the ansatz 

f a[0] = / d^'x 



I {dct^f + \m'cf' - Va(p) 



where (we drop the subscript in the following) V{p) is the effective potential and p = \ 
a field independent term, the A-dependence of V may be derived from the flow-equation 
one finds, using for the cutoff function 9(|fc|, A) = 6(|A;| — A), ||T5| 



A 



dV{p) 
dA 



^ A3 , 
-T— In 

27r2 



1 — exp 




e(A2 + V' + 2pV" 



(4) 



(5) 

Up to 
m and 



(6) 



This equation has to be solved numerically. One way to proceed is to expand the effective po- 
tential around some field-value p, truncate the expansion at some finite order, and thus obtain a 
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finite system of coupled ODEs for tfie "couplings" V^'^\p) (tliis approacli is known as tlie "local 
polynomial approximation" or LPA). This approach was followed for the present flow-equation in 
||15|| , where values for the critical exponents 6, u, and 1] where obtained (see table below). In the 
ERG-formulation in three dimensions, critical exponents where also derived along these lines 0. 
However, it turns out that the results obtained in the TRG-approach using the LPA are relatively 
poor compared to the results obtained from lattice studies or higher order calculations in the 
e-expansion. This should not be connected to the derivative expansion, where we expect errors of 
the order of the critical exponent t], being ^ 0.03 for the one-component theory. 

In this letter we will thus study the numerical solution to the full partial differential equation 
(^. This has the further merit that we will be able to give the approximate equation of state, 
namely the function "^^Jf^ = (f)V'{p) in the vicinity of the critical temperature. We use the 



algorithm proposed in ||30|: We discretize in p and use the flow-equations for V\p) and V"{p) 
at each point. The flow-equations involve the derivatives of V up to V"*^^-*, and we obtain these 
derivatives from matching relations by imposing continuity of V and V" at the points. In the 
(low temperature-) broken phase, it turns out convenient to also follow the A-dependence of the 
minimum of the potential, given by 



A 



dA 



47r2 



3 + 2 



(3)~ 



A2 + 2poV^' exp V'^ + ^ 



1 



-e(A2 + 2poK') 



(7) 



Furthermore, for a study of the critical properties at a second order phase-transition, one should 
make use of the fact that properly rescaled quantities asymptotically exhibit scaling at the critical 
temperature. We introduce the three-dimensionally rescaled potential and field through 



P 

AT 



u 



V 



d'V 
dp^ 



Using A = ^, the flow-equations for u and now read 



du 
^dX 



-3u + u'k In (l 

27r^ 



exp 



and 



dn^ 



(8) 



AVl + u' + 2ku"] } 6(1 +u' + 2ku") 



(9) 



1 



1 + 2kqu'^ 



exp XJl + 2kquU - 1 



Q{1 + 2kouI) (10) 



The three-dimensional rescaling in (|]) is motivated by the fact that for A^ + V' + 2pV" <^ T^, 
i.e. when both the external scale and the mass-scale of the theory are small compared to the 
temperature, the flow-equation for m, eq. (^, reduces to the flow-equation of the purely three- 
dimensional theory with a 0-function cutoff (up to a field- independent contribution): We have 



du 
^dX" 



-3m- 



-u'k- 



1 



4vr2 



\n{l+u' + 2ku") 



27r2 



In A 



Q{l+u' +2k,u") + O 



^A^ + V'+2pV" 
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Figure 1: The dimensionless minimum of the effective potential as a function of A/T for different 
temperatures above and below the critical one ((7(Ao) = 0.1). The solid line gives the scaling 
solution of the purely three-dimensional theory. Dimensional reduction is nicely observed. 



This is the manifestation of dimensional reduction and the dominance of classical over quantum 



fluctuations in the critical domain [pl[|. 



3. We now proceed to the results for the critical behavior of the theory obtained by a study of 
the system of equations resulting from (||) - ([T0|). As was already pointed out in ||I5|, the phase- 
transition is found to be second order, consistent with the expectations from other theoretical 
studies and experiments. We will in this letter concentrate on the critical exponents and the 
equation of state. More detailed results on other aspects of the critical behavior like critical 



amplitudes and nonuniversal quantities will be presented elsewhere |^2[- a starting potential 
at Ao ^ T (in practice we use Aq = exp(3.5)T ^ SOT) one needs the effective potential of the 
theory at vanishing temperature. We use here the tree-level potential 

^Ao(p) = ^(p-Po(Ao))' (12) 

and a small value of 5'(Ao) = 0.1. In principle one could use the one- or two-loop potential, but 
this would only change the nonuniversal results. Below, all dimensionful quantities are given in 
units of y/po which we arbitrarily set to 1. 

At the critical temperature, the theory exhibits scaling: The dimensionless minimum of the 
potential approaches a nonvanishing constant as A = ^ ^ 0. Figure 1 displays Kq as a function 
of A for different temperatures above and below the critical temperature. Above the critical 
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This worli 
TRG+LPA 


15| 


0.35 


1.32 


5.0 
3.57 




0.015 


0.66 (= Z//2) 
0.58 


3d ERG + LPA § 
3d ERG [|] 


0.333 
0.336 


1.247 
1.258 


4.75 


0.045 
0.044 


0.638 
0.643 


Best values 




0.325 


1.240 


4.81 


0.032 


0.630 



Table 1: Critical exponents from different approaches 



temperature the theory is in the symmetric phase and the minimum reaches for finite values of 
A. Below the critical temperature, the theory is in the broken phase and the dmenszon/u/ minimum 
Po = AT^Ko reaches a constant as A — > - the dimensionless minimum diverges. Also included 
is the result of a study of the effective three-dimensional theory using the leading contribution of 
([TT]) only. In this case, the parameter controlling the deviation from the critical line is not the 
temperature but the value of the minimum at the starting scale. One nicely observes the effective 
three-dimensional behavior of the finite-temperature theory close to the phase transition after an 
initial running which is governed by the four dimensional theory (in the case of the TRG the 
dimensionful parameters are constant in this regime, since the quantum fluctuations are included 
in the starting action already) for A 3> T and a "crossover" regime for A «i T. The TRG in this 
way provides the possibility of studying dimensional reduction nonperturbatively. 

The existence of the asymptotic scaling solution implies universal critical exponents and a 
scaling equation of state. From the study performed in the present work we can extract all critical 
quantities directly in the full finite-temperature theory. In the following we concentrate on the 
exponents /3, 7 and b and the equation of state. The approximation performed by neglecting 
nontrivial terms in the derivative expansion prevents us from studying the exponent 77 which is 
given by the anomalous dimension of the order-parameter field at the critical point. We thus 
by assumption have = 0, which has to be compared to the best value in the literature, giving 



r\ = 0.032 p5|. As stated above, the approximation is rather modest in the case under study. We 
expect the other critical exponents to have an error of the order of 77, i.e. of about 5 %. 

The exponents we will consider are defined in the following way: The exponent /3 describes 
the temperature-dependence of the magnetization (in our case of the expectation value of 79) as 
the temperature approaches the critical temperature from below: (yjo = v^2kJ^) 

¥Pocx(r,-Tf (13) 

The exponent 7 describes the behavior of the susceptibility, in our case corresponding to the 
inverse mass of the order-parameter field, as T reaches T^. We have 

V'(po)+2po^"(po)(x|Te-Tr (14) 

Finally, the exponent 8 describes the form of the effective potential for small p ai T = Tc through 

V'ip) oc p^ (15) 
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Figure 2: V\p) in the critical regime from the thermal renormalization group (diamonds). For 
comparison we give the properly rescaled results of a lattice study and the e-expansion to 
order e°, and |]33 |. 



Table 1 gives our values for the critical exponents in comparison with results from other studies. 
There is a clear improvement of the results compared to the exponents found using the LPA in 
T5[l . Compared to the best known values, the errors are of the expected order. The results could 



be further improved by going beyond leading order in the derivative expansion. 

We finally turn to the equation of state. This has been addressed in the three-dimensional 
theory through the e-expansion, ERG-calculations, and lattice studies. We will compare our results 
below. The equation of state is essentially given by the derivative of the effective potential as a 
function of the order-parameter field (in a solid-state context the magnetic field if as a function 
of the magnetization M) in the critical regime. Up to the normalization of the potential and the 
field, the equation of state is universal, and may for not too large values of the field be computed 
in the three-dimensional theory (for very large field-values, dimensional reduction generally does 
not hold since the mass of the order-parameter field is large compared to the temperature). We 
choose here a temperature slightly below the critical one but safe inside the critical domain and 
plot V\p) in figure 2 (diamonds). For comparison, we also display the properly rescaled results 
of a lattice study , which is up to normalization given by the simple expression 

^ ' -1 -4p + 12p2^ 



y\p) 



(16) 



as the solid line as well as the results found from the e-expansion to order e°, e and e^ |3[ (the 
dashed lines labeled by (0), (1) and (2) resp.). This figure may be compared to figure 5 in p6| , 
where also a comparison of the lattice results and the results of the ERG-study in the three- 
dimensional theory [H is given. One observes an impressive agreement of both the ERG-results of 
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1^ and the TRG-results presented here with the lattice-calculations, while the e-expansion gives 
accurate values for the critical exponents but works only rather modestly for the equation of state. 

4. In conclusion, we have presented an approximate solution of the thermal renormalization 
group-equation for the one-component selfinteracting scalar theory in the critical domain. We 
work in lowest order in the derivative expansion but make no approximations to the form of the 
effective potential. We explicitly obtain a nonperturbative description of dimensional reduction 
and demonstrate how the critical behavior of the theory is dominated by three-dimensional classical 
fluctuations (figure 1). The critical exponents corresponding to the second order phase-transition 
that the theory exhibits are obtained with an accuracy of the order of 5 % already within the 
relatively simple approximations performed here and could be improved by going to higher orders 
in the derivative expansion (table 1). We finally find results for the equation of state of the theory 
that are in surprising agreement with the results of lattice studies (figure 2). 

All in all this work demonstrates that the thermal renormalization group is a valid tool for the 
study of critical phenomena at finite temperature. It has a number of advantages compared to 
other methods, maybe the most prominent one being the fact that it is manifestly gauge- invariant. 
It thus simplifies the study of the critical behavior of gauge theories as compared to the usual 
ERG-approach since no non-gauge-invariant counter-terms are enforced through the introduction 
of the cutoff. Furthermore the TRG makes it possible to work with effective actions for gauge- 
theories also for nonvanishing values of the external scale, corresponding to the course-grained free 
energies used in statistical mechanics for the study of first order phase transitions. Applications 
to the critical behavior of gauge theories should thus be very interesting. 

Acknowledgment: I would like to thank Jiirgen Reingruber, Chris Stephens and Johannes 
Manus for interesting discussions. 
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